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Abstract 

Motivated by the search for gravitational waves emitted by binary black holes, we investigate the 
gravitational radiation field of point particles with spins within the framework of the multipolar- 
post-Newtonian wave generation formalism. We compute: (i) the spin-orbit (SO) coupling effects 
in the binary's mass and current quadrupole moments one post-Newtonian (1PN) order beyond the 
dominant effect, (ii) the SO contributions in the gravitational-wave energy flux and (iii) the secular 
evolution of the binary's orbital phase up to 2.5PN order. Crucial ingredients for obtaining the 
2.5PN contribution in the orbital phase are the binary's energy and the spin precession equations, 
derived in paper I of this series. These results provide more accurate gravitational-wave templates 
to be used in the data analysis of rapidly rotating Kerr-type black-hole binaries with the ground- 
based detectors LIGO, Virgo, GEO 600 and TAMA300, and the space-based detector LISA. 

PACS numbers: 04.30.-w, 04.25.-g 
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I. INTRODUCTION 



The aim of this paper is to derive the spin-orbit coupling terms in the gravitational radia- 
tion field of compact binary systems one post-Newtonian (1PN) 1 order beyond the dominant 
effect. Paper I of this series [l| dealt with the problem of the spin-orbit contributions in the 
compact binary equations of motion at 1PN relative order. 

Our motivation is the on going search for gravitational waves (GWs) emitted by in- 
spiralling binary systems of spinning, and possibly maximally spinning, black holes in the 
network of detectors LIGO (Laser Interferometer Gravitational Wave Observatory), Virgo, 
GEO 600 and TAMA300, and the future search with the space-based detector LISA. When 
the Kerr black holes are maximally spinning (or close to maximal), the GW templates need 
to take into account the effects of spins, not only for an accurate parameter estimation 0- 



but also for a successful detection [8Hl9|. Furthermore, spin effects should be included at 
PN orders beyond the currently known dominant spin-orbit and spin-spin terms. The con- 
tributions of spins are added to the templates developed for the case of non-spinning binary 



black holes or neutron stars, and which are currently known at 3.5PN order [20H23|. The 
spins represent some of the possible effects depending on the internal structure of the bodies 
which can be numerically important in the LIGO /Virgo bandwidth. This is true even if we 
observe only the inspiral phase of moderate-mass black holes with individual mass less than 
10 M Q . Within the PN formalism the compact objects are treated as point particles. It is 
then natural to model spinning black holes as point particles with spins. 

The equations of motion including the spin-orbit (SO) effect were obtained in paper I at 
1PN relative order, which corresponds formally to the 2.5PN order beyond the Newtonian 
force law in the case of maximally rotating compact objects. Paper I essentially confirmed 
the equations of motion derived previously by Tagoshi, Ohashi and Owen [24|. Furthermore, 
paper I derived the complete set of Noetherian conserved integrals of motion at that order 
(namely, 2.5PN for the spins). In the present paper, we tackle the problem of the gravi- 
tational radiation field at the same 2.5PN order, using the multipolar PN wave generation 



formalism of Refs. [25H30j . More precisely, we shall compute here the SO contributions in 
the compact binary's mass-type and current-type quadrupole moments, both of them with 
1PN relative accuracy. These moments, together with some easily computed higher multi- 
pole moments which necessitate only the lowest-order precision, are necessary to compute 
the total GW energy flux J 7 . The computation of the current quadrupole moment was 



previously attempted in Ref. [31], but we shall point out two important flaws in that refer- 
ence (see below for details). Our result for the current quadrupole moment is substantially 
different from that of jHf. Concerning the mass quadrupole moment, it is computed here 
for the first time. Having in hand the total energy flux, using the center-of-mass energy 
E computed in paper I, we deduce (by energy balance arguments) the equation of secular 
evolution of the binary's orbital frequency. The latter is the crucial ingredient needed to 
build GW templates for spinning compact binaries. 



To describe particles with spins we use the formalism originally developed in Ref. [32H35 



(an effective field theory scheme has recently been proposed [36-38]). This formalism has 



already been successfully applied to the problem of spinning compact binaries in Refs. [24, 



, and, in the test-mass limit case, in Refs. HEl. In particular, Kidder, Will and 



1 By nPN we refer to the terms of relative order (v/c) 2n where v is the binary's orbital velocity and c the 
speed of light. 
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Wiseman derived in Ref. 39|, |40| the lowest-order spin-coupling effects — at 1.5PN order in 



the case of maximal Kerr black holes — , and the first spin-spin effect, quadratic in the spins 
- appearing at 2PN order — , in the equations of motion and the gravitational radiation 
field. As we shall see below we find complete agreement at that order with their results. 



The present paper together with paper I extend therefore the works [39|, |40[ to include the 
next-order spin effects. Since those effects, of 2.5PN order, are linear in the spins (the next- 
order spin-spin term coming along at 3PN order), we complete the derivation of all the spin 
contributions in the GW form up to 2.5PN order. 

The paper is organized as follows. In Sec. [XT] we review the general formalism for wave 
generation from arbitrary PN sources. In Sec. II I II we compute the multipole moments 
of compact binary systems at the lowest PN level in the spins (which means 1.5PN for 
mass moments, and 0.5PN for current ones). Sec. IIVI constitutes the core of the paper. 
We compute there the mass and current quadrupole moments at 1PN relative order, i.e. 
2.5PN and 1.5PN for the mass and current types, respectively. In particular, the crucial 
contribution of non-compact-support terms (which are sourced by the gravitational field 
itself) is obtained. The final results for the multipole moments and the GW flux are presented 
in Sec. [V] and in Sec. |VT1 we reduce those results to the physically relevant case of quasi- 
circular orbits. In Sec. IVHt we express the main equations defining the binary evolution and 
the GW signal in terms of spin variables with constant magnitude, which generalize those 



used in Refs. 39|, |40). In Sec. IVIIIt relying notably on the results of paper I, we obtain 
the secular evolution of the orbital frequency in the case of circular orbits and discuss some 
implications for ground-based and space-based detectors. Section [IX] summarizes our main 
conclusions. 

All the notations and conventions are the same as in paper I. Notably, in order to explicitly 
display the powers of 1/c appropriate for maximally spinning compact objects — in which 
case the spin variable is formally of the order 0.5PN ~ 1/c — , it is convenient to adopt 
as basic spin variable a quantity having the dimension of an angular momentum multiplied 
by c, namely S = c S tTue (see paper I for discussion). The precise definition of the spin 
contravariant vector we use and the supplementary condition it satisfies are given in Sec. II 
of paper I. See the same reference for the details concerning the stress-energy tensor of 
spinning point particles. 



II. POST-NEWTONIAN SOURCE MULTIPOLE MOMENTS 

We start with a brief review of the PN multipole moment formalism at the basis of this 



approach (full details can be found in Refs. [2l|, |28|-|30j). This formalism is valid for general 



localized material sources, satisfying the usual PN requirements of weak self-gravity, slow 
motion and weak internal stresses. In particular, the size of the source a has to be small 
with respect to the typical (reduced) wavelength X of the gravitational radiation this source 
produces, i.e. a/X = 0(e), with e ~ v/c being the slowness PN parameter. We shall 
abbreviate it as e = 1/c and denote the PN remainder terms by 0(l/c n ) henceforth. 

Let = (ct, x) be an harmonic coordinate system covering the whole material source. 
We pose = \f—g g^ w — r] fiu where g^ v and g are the inverse and the determinant of the 
usual covariant metric g^ u , and where if"" denotes an auxiliary Minkowskian background 
metric (Greek letters represent space-time indices; our signature is +2). The Einstein field 
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equations, relaxed by the harmonic-coordinate condition, d u h fMU = 0, take the form 

nhflv = lfwrG^ _ M | p | T ^ + A ^^] ? (2.1) 

where □ = ri pu d^ v is the flat d'Alembertian, and where the second equality defines the 
stress-energy pseudo tensor r pu of the matter and gravitational fields in harmonic coor- 
dinates. Here T pu is the stress-energy tensor of the matter fields and K pu [h] represents 
the gravitational source term, namely a complicated non-linear functional of h pa and its 
space-time derivatives d\h pa (see e.g. j2l[ for the expression). We shall see that K pv gives a 
crucial contribution to the multipole moments at the relative 1PN order in the spins. The 
stress-energy pseudo tensor is conserved by virtue of the harmonic-coordinate condition, 

d u h^ = => d u T^ = 0. (2.2) 

The multipole moments of the source are generated by the components of the pseudo 
tensor T pu or, more precisely, of its formal PN expansion r pu = PN[t^]. In this sense, 
the formalism is physically valid for PN sources only. The PN expansion has a special 
structure which can be matched to the exterior multipolar field of a PN source, allowing one 
to define an appropriate notion of PN multipole moments It is convenient to define 

(Latin letters representing space indices) 

£ = C - 2 [r 00 + t 11 ] where T u = %7" , (2.3a) 
E^eec- 1 ^, (2.3b) 
E {j = r ij . (2.3c) 

The mass- type moments Ii{t) and current-type ones Jz,(i) are referred to as the source 
multipole moments in order to distinguish them from the so-called radiative moments, seen 
at infinity and generally denoted Ui{t) and V/^t). They are given by 2 

h (t) = FP / ri 3 x |x|» /' iz { W i L S - c2( ^[^ + 3) W«) ** % 



2(2^ + 1] 



c\i + 1)(£ + 2)(2£ + 5) 
J L {t) = FP s ab (i e J rf 3 x|x| B y dzls e (z)x L - 1 ) a E b 



Se+2(z) (x,t + z\x\/c) , (2.4a) 



2£ + 1 



c 2 (£ + 2)(2£ + 3) 



S i+1 (z) Xl-i)oc S 6c |(x, t + z\x\/c) . (2.4b) 



These expressions are general, in the sense that they are formally valid at any PN order, so 
that there are no remainder terms 0(c~ n ) involved. The dots indicate the time derivatives 



2 In our notation, L = i\ ■ ■ ■ ig represents a multi-index composed of I multipolar indices - ■ ■ ,it, and 
Xl = Xi t ■ ■ ■ Xi e stands for the product of £ spatial vectors x l = Xi. We denote the symmetric-trace- free 
(STF) projection by means of a hat over the tensor symbol, xl = STF(xi ± •••Xi t ), or of brackets () 
surrounding the indices, xl = X (L)- When the multi-indices L are summed up (dummy indices), we omit 
to write the I summation symbols from 1 to 3 over their indices. 
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and E a bc is the Levi-Civita symbol in 3 dimensions with £123 = 1. Notice the peculiar feature 
that, besides the usual spatial integration, the moments involve an extra integral over a 
variable z defining a "cone" of integration u = t + z|x|/c; hence the sources depend 
on the point (x, u) as indicated. The "weighting" function associated with the ^-integral, 
8e(z) = ^t+n\ (1 — z 2 ) e , is normalized in such a way that J ^ dz5e(z) = 1. When performing 
explicitly the PN expansion of the moments, the ^-integration is to be transformed into an 
infinite local series (in the sense that z is integrated out), which constitutes the basis of the 
practical evaluation of the multipole moments. Namely, we have 



1 +2? (of + IV /Ixl d\ 2k 



~ 1 k=0 



The above expression is then inserted into the right-hand side (RHS) of Eqs. (12 .4p . and 
truncated to fit with the PN order of the calculation. 

A crucial finite part (FP) procedure is involved in the definition of the source multipole 
moments (12.41) . It consists of (i) multiplying the integrand of the moments by a regularization 
factor |x| B where B is a complex number, 3 (ii) performing the Laurent expansion when B 
tends to the "physical" value B = 0, and (iii) picking up its finite part FP, namely the 
coefficient of the zeroth power of B. The finite part regularization is therefore equivalent 
to removing the poles B~ x , B~ 2 , ... (in the analytically continued -B-dependent integral) 
before taking the limit B — > 0. The FP procedure is needed to compute the non-linear 
contributions to the moments (generated by the gravitational source term A My ), which have 
a non-compact support extending up to spatial infinity. Notice that no assumption nor 
physical restriction (in principle) is involved, in the latter FP procedure, which has been 
proved [28-3(j to yield the correct expression of the multipole moments for general extended 



PN sources. It is precisely the FP that guaranties this within the present formalism; the 
divergent terms B^ 1 , B~ 2 , ... have no direct physical significance. Such FP when B — > is to 
be carefully distinguished from the self- field regularization (e.g., Hadamard's or dimensional 
regularization) which is to be invoked when treating the application of the general formalism 
to singular point-particle sources. 

We emphasize that the expressions (12.41) constitute a priori only a definition of the 
source multipole moments. The point is that such definition is fully related to the physical 
asymptotic wave form at infinity from the source, which is computed using the multipolar 
post- Minkowskian formalism 25, 26[. In particular, the radiative multipole moments U^t) 



and Vi,(t) which parameterize the asymptotic wave form are given by some non-linear func- 
tionals of the source moments (12. 4p . made of many non-linear interactions between them, 
including the famous GW tails corresponding to the interaction between these moments and 
the total monopole mass M of the source. The tails have been computed within this ap- 



proach in Refs. [46|, |47|. However, for deriving the spin terms at 1PN relative order, all these 
non-linear multipole interactions, notably all the tails, are negligible (see further discussion 
below). It is therefore sufficient to consider only the source multipole moments of Eqs. (12.41) . 



3 Generally the regularization factor is taken to be |x/ro| B where rg denotes some arbitrary constant length 
scale. Here we set ro = 1 for convenience. 
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III. LOWEST-ORDER SPIN EFFECTS IN THE MULTIPOLE MOMENTS 



The multipole moments discussed above can be applied to any source. Here, we specialize 
them to binary systems of point-particles with spins. The formalism to describe particles 

and constitutes the basis of most subsequent 
We reviewed this formalism in paper I and 



39 



with spins was developed in Ref s. J32U35 
computations in this field y/ 
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refer the reader to this paper for details and notation. The spin contribution (marked by 
the underneath label S) to the stress-energy tensor of the particles reads 



T^(t,x) 

s 



-V 



<Hx-yQ 



1 -H- 2 



(3.1) 



where 5 is the Dirac three-dimensional delta-function, V p denotes the covariant derivative, 
v i{t) = ( c i v i) with v\ = dyl/dt being the coordinate velocity of particle 1, g± is the deter- 
minant of the metric evaluated at the location of particle 1 (following Hadamard's self-field 
regularization), and 1 -H- 2 means the same expression as preceding, but for particle 2. The 
anti-symmetric spin tensor S± u (t) is introduced in Sec. II of paper I. The covariant four- 

^ £ nupa u i<gu ^ j g transverse to the particle's four- velocity, 



vector S,, is defined by S, — — _ "„■-■„■ 
S^Ui = 0. All results below are expressed in terms of some particular space-like contravari- 
ant spin variables for each of the particles, namely S{ and SI, the definition of which can be 
found in Eq. (2.19) of paper I. 

Similarly to the quantities introduced in Eq. ( 12. 3p . we define the following matter- 
source densities, depending on the components of the spin stress-energy tensor (13. 1ft : 



cr = c 

s 



rpOO _|_ rpi 



with V" d u •/•" . 



Oi = c- 1 T 01 , 
s s 

°Ai = T ij • 

s s 



(3.2a) 
(3.2b) 
(3.2c) 



They are such that their "non-spin" counterparts (say ns admit a finite non-zero limit 
when c — > +oo. They read: 



a 



o = -^e tjk vlS{d k 6 1 + 1^2 + o(\ 

S C y C 

<Ji = -^EijkSidA + 1 -H- 2 + , . 



(3.3a) 
(3.3b) 
(3.3c) 



where we denote Si = 5(x — yi) and where d k 5i means the gradient of 5\ with respect to 



the field point x 



.r 



). As shown by Eqs. (I3.3p . the leading order of the vector and tensor 



densities sO"i and s^ij is 0.5PN ~ 1/c. However, the scalar density s<r starts at a higher 
level, being at least 1.5PN ~ 1/c 3 . At leading order in the spins, the E^'s, which depend 
on the contributions of both matter and gravitational fields according to Eqs. (12. ip and ( 12. 3p . 
will reduce to their compact-support material parts, namely the s^i/'s given by Eqs. ( 13. 3p . 
Indeed, the non-compact support gravitational part, whose origin lies in the source term 
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K^ u present in the RHS of the field equations (12.11) . always appears at a sub-dominant level, 
1/c 2 beyond the leading PN order. Hence, 

S-a+ofiV (3.4a) 



s s v c 5 



Si = & + O ( \ ) , (3.4b) 



s " s • V c 3 



Ey = a y + O (±) . (3.4c) 

The non-compact-support gravitational source terms play a role in our computations at the 
next-to- leading order only (see Sec. II Vj) . We conclude that the dominant contribution to 
the multipole moments (I2.4p due to the spins is given by 

/, = / {fa g - ^"ff + 3) «* f .} + O (1) , (3.5a) 

Jl = £o6(ii y rf 3 xx L _i >a + C (^A . (3.5b) 

Thus, the dominant order is 1.5PN ~ 1/c 3 for spins in the mass-type moments II, but 
only 0.5PN ~ 1/c in the current-type ones Jl- It is then evident (mathematically and 
physically) that the spin part of the current moments always dominate over that of the mass 
moments. We insert the explicit values (I3.3P into Eqs. ( I3.5p . integrate in a straightforward 
way (resorting to an integration by parts and using the properties of the delta-function) and 
get 

tv\ S{ e im -(£-1) y{ S{ e m y^ 2) 



l L c 3 (£+l 

+ 1^2 + OK] . (3.6a.) 



Jl = -7T- yr 1 S\ e} + 1^2 + - . (3.6b) 



_ ±^yt l S^ + l^2 + 0{\ 
s 2c \c 3 

It is worth to mention that in this calculation, limited to the lowest PN order, the spins can 
be considered as constant since their time variations, as given by the precessional equations 
(see paper I), are always smaller by a factor 1/c 2 at least. 



IV. HIGHER-ORDER SPIN EFFECTS IN THE MULTIPOLE MOMENTS 

For the present purpose, we need the spin contributions to the mass-quadrupole moment 
Jjj and the current-quadrupole moment one PN order beyond the leading terms obtained 
in Eqs. (13. 6p . This means at 2.5PN order for 1^ and at 1.5PN order for J^. As said pre- 
viously, the non-linear gravitational source terms, with non-compact support, start playing 
a role at the 1PN relative order. 4 Therefore, they do make a net contribution to the spin 



4 The rule admits some exceptions, though. For instance, in the non-spinning part of the mass multipole 
moments II, one may have expected the non- linear non-compact gravitational terms to appear at 1PN 
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parts of both at 2.5PN order and Jjj at 1.5PN order. We can note here that the authors 
of Ref. 31] computed at 1.5PN order but neglected all the non-compact source terms in 
their calculation, thus obtaining an incorrect result. 

We now reduce the multipole moments to the required PN order by inserting the expan- 
sion formula (I2.5P into the general expressions (I2.4p . Neglecting PN terms of higher order, 
as indicated, we have 



FP 

B=0 



cTx |x 



B 



1 A 
14^' |X| 



504c 4 ^ 



' 2i c 2 Xi i k ^ k 



10 



L89c 4 " 
B 



%ijk | x | 2 Sit + ——, r Xijk&ki \ + O ( — 



= f£ ^ / d X ' X ' i Sb + 14c2 £j> 



54c 4 ^ 

\ 2 t h - 



28c 2 



j')ac 2-*bc 



+ o 



(4.1a) 
(4.1b) 



where FP denotes the essential process of extracting the finite part, as explained in Sec. [TTJ 
We can further reduce Eqs. (14. ip by substituting the appropriate expressions of the source 
terms £ Mi , as given by Eqs. (4.11) of Ref. |28J. For completeness, we list all the necessary 
expressions below: 



4V 2 - 



a 



ttGc 2 



diVdiV 



+ 



i r 



-Vd 2 t V - 2Vid t diV - WijdijV - -{d t Vf 



7rGc 4 \ 

+ IdiVjdjVi - diVdiW 
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VdiVdiV 



E 



1 + 



AV 



Gc 2 



7T 



d k V(d t V k - d k Vi) + -dtVdiV } + O 



4 



(Tin H 

y 4ttG 



d i Vd j V--6 ij d k Vd k V\ + 



(4.2a) 
(4.2b) 
(4.2c) 



where the material source densities are given by a = c [T + T M ] , Cj = c T and 
crjj = T^. The non-compact support terms in Eqs. (14.21) are parameterized by a particular 
set of "elementary" potentials V, Vi, Wij, which enter the harmonic-coordinate near zone 
metric at the 2PN order computed in Ref. 48]. 5 Their complete expressions are given in 
Sec. Ill of paper I. 

We can simplify S substantially by using some identities of the type diAdiB = 
| [A (AB) — AAB — BAA] , replacing the Laplacians AA and AB by their PN sources, 
and disregarding the pure Laplacian term |A (AB) because it makes zero contribution to 
the moment. This last point comes essentially from the fact that, after integration by parts, 
a pure Laplacian term in the moment sources is equivalent to a source term proportional to 



order, but these terms turn out to be in the form of a pure divergence and can be integrated out to zero. 
As a result, the non-compact terms of the non-spinning parts contribute to II at 2PN order only (and at 
1PN order to J L ). See Refs. [21 H for details. 

The potential called Wij in Ref. 28] differs from the present Wij (whose definition is given in paper I) 



according to the formula Wij = Wij — \8 l jW ' , hence W = Wu = — \W 
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some Axl = 0. Beware that, because of the presence of the finite part FP, such a "proof" 
is not correct in general and, in fact, the pure Laplacian terms do generally contribute at 
high PN orders. However, for the terms under concern, merely at the 1PN relative order, 
the argument can be made rigorous and shown to work, so that we can indeed discard these 
pure Laplacians in the present computation (see Ref. H| for the proof). Hence, we get the 
simpler formula, 

S ' = " + + [-2V,dAV - WijdijV - l -(d t V) 2 + TdiVfijvX 

+ O (1) , (4.3) 

which differs from E by pure Laplacian terms of the type ~ A (AB) . We have checked that 
the two different forms £ and Eqs. (14. 2 aft and ( 14. 3D . lead to the same final result. 



A. Compact support contribution 

Having now the general set up for our computation we consider first the compact-support 
part of the multipole moments, i.e. that part proportional to the material source densities 
a^ u , and given by the first terms in Eqs. (I4.2p - (l4.3p . For these terms we make two compu- 
tations. The first one consists of (i) evaluating all the components of T^ v to the correct PN 
order [extending thus Eqs. ( 13. 3ft and including the monopolar sources, which may give rise 
to spin terms through the metric], (ii) computing their time derivatives by making use of the 
usual replacement of accelerations by the equations of motion, and of the time-derivatives 
of the spins by the precessional equations, (iii) transforming the time-derivatives, when ap- 
plied to delta-functions, into spatial ones using the formula d t 8\ = —v\di5i, (iv) operating 
by parts the spatial derivatives of delta-functions to finally integrate thanks to the basic 
property of delta-functions. 

Normally, such basic property of delta-functions reads J rf 3 xF(x)5i(x) = (F)i, where 
(F)i is simply the value of the function at the point 1. When the function is regular, there 
is no problem and we have (F)i = F(yi,t), for instance (xl)% = y\. However, when the 
function F is singular at the point 1 {i.e. when x — > yi), a choice must be made for 
a "self-field" regularization, able to subtract the infinities in a consistent way. There are 
various possibilities. In the present work we adopt, following paper I, the Hadamard self-field 
regularization. At the order we are working (relative 1PN order) the various possible choices 
are equivalent. For instance, one can show that dimensional regularization would give the 
same result as Hadamard's regularization, essentially because at such low PN order there 
are no poles in the dimension of space [say, oc (d — 3) _1 ], which correspond to logarithmic 
divergences in Hadamard's regularization. We then define (-F)i to be given by the partie 
finie of the function F at point 1 in the sense of Hadamard (see e.g. [49| for a full account 
of this regularization). Suppose for example that F = U xl where U = ^r 1 - + is the 
Newtonian potential of point particles, singular at the location of the particles. Then we 
easily compute that Hadamard's partie finie is {F)\ = ^f-yf- 

Our alternative computation of the compact-support terms in the moments is the same as 



the one performed by Owen et al. 3J|. It consists of applying the following formula (derived 
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Ill 



t) T» u (t, 
s 



d 

cdt 



0(M v) 



+ 



1 per °1 



J- op D i ^l 



1 2 



(4.4) 



which is valid for any function F(x, t), and where (F)i and (d p F)i have to be understood 
as the Hadamard partie finie of the function and its derivative. This formula is very useful 
but must be handled with care. In particular, when computing (d p F)i in the second term 
of the RHS of (14.40 . we notice that the gradient is to be taken first, and only then should 
one deduce the value at point 1. The result can be different if one permutes the order of 
operations. Suppose for instance that one is computing (d t F)i where F = xl. Clearly, since 
d t XL = the result is zero. However, if one computes 4x{F)\ instead of (d t F) 1 , one obtains 

an incorrect non-zero result, which is equal in this case to ^ \Vx\ = ty[ 



We found 



that this error, 



i.e. computing %{xl)i = £y[ L 1 v± instead of {d t XL)i = 0, was committed 
in the evaluation of the (compact-support part of the) current quadrupole moment Jjj in 
Ref. (E . 



B. Non-compact support contribution 

We now derive the non-compact support part of the multipole moments. Inspection of 
the expressions (I4.2p -f l4~3]) shows that we need only the elementary potentials V , Vi, 
(and W = Wa) at their lowest PN order in the spins. They read: 

V = ~ e ijk v\ S{ d k + 1 o 2 + O f^j , (4.5a) 
Vi = ~ e ijk S( d k (±) + 1 ^ 2 + O (i) , (4.5b) 
Wij = -- c end v{ ] S? d l (^j + | <% e klm vi S{d m U-J+1^2 + f^j , (4.5c) 



W = —e klm v k S l 1 d m [-)+1^2 + 0[-) . (4.5d) 



— e klm v k 1 S[d m (-)+1^2 + o(- 
s c \ r i/ V c ' 

The latter spin contributions enter only the "non-spin" parts of Eqs. (14.20 — (14.30 . We need 
now to perform the integrations in Eq. (14.11) . The calculation becomes simple if we use 
the following trick. In the multipole moment's integrands we transform all the gradients 
evaluated at the field point x into gradients evaluated at the source points yi or y 2 using 
e -9- -£j (l/ r i) = — ~§^f (V r i)- Then, we put the source-type gradients and outside 
the integrals and we express the result solely in terms of the function 

Y L ( Yl , y 2 ) = - J- FT f d 3 x |x| B ^- , (4.6a) 



27T b=o J r x r 2 
which is known to admit the analytically closed-form [2 



Yr 



^E^"^. (4.6b) 
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Thus, the closed-form expressions of the non-compact (NC) parts of the spin multipole 
moments (they depend on the function Yj, for £ = 2,3) are: 



(NC) 

s ij 



2Gmo 



£mnk ^1 Q n dkl^ij ^mnp ^1 Op A Yi 



1 2 



1 2 



^kmn^nSi dndkl^ij 2£i mn Do Si dkndlYij 
12 12 



+ 21 lmn 1 dt 



1 



<9nfc dlYijk 



jjj rim 

21 fcmn 1 dt 



<9m" 



1 2 



+1«2+0|- 



(NC) 



Gm 2 



^kmn Si dk dlnYj)k 4" ^Imn <Jl <9p dpnXj)k 
2 1 ' 2 1 ' 



+ 1 <R> 2 + C — 



(4.7a) 



(4.7b) 



where <9ifc = d/dyf and <9 2 fc = djdy\. The final computation of the moments using for- 
mula (I4.6bl) is straightforward. 



V. RESULTS FOR THE MULTIPOLE MOMENTS AND FLUX 



The expressions for the multipole moments $1^ and sJiji including both compact and 
non-compact contributions as computed in Sec. IIVI are quite long if written in a general 
frame. They can be substantially simplified by going to the frame of the center-of-mass 
(CM). When working in the CM frame it is convenient to use the following spin variables: 



S 

£ = m 



Si + s 2 , 
S 2 _ Si 

m 2 mi 



(5.1a) 
(5.1b) 



These spin variables were initially introduced by Kidder |40| except that here we denote 
£ what he calls A. Mass parameters will be denoted by m = mi + m 2 , 5m = mi — m 2 
and v = mi m 2 /m 2 for a mass ratio such that v — 1/4 for equal masses and v — > in the 
test-mass limit. 

The CM frame is defined by the nullity of the binary's dipole moment or equivalently the 
CM vector G. At 2.5PN order including spin effects, it can easily be determined using the 
vector G evaluated in paper I. However, here we need only the lowest order term (1.5PN in 
the spins) together with the 1PN non-spin correction; the 2.5PN term in the spins cancels 
out. To the needed order we have (see e.g. Eq. (5.13) in 



Yi 

y-2 



m 2 v 5m 

m 2c 2 m 

mi v 5m 

m 2c 2 m 



v 2 - 



Gm 
r 

Gm 

r 



x + 



V 



m cr 



mc 



v x S. 



3 vxS, 



(5.2a) 
(5.2b) 



which gives the CM positions of the particles, yi and y 2 , in terms of the relative position 
and velocity, x = y x — y 2 and v = dx/dt = vi — v 2 (we pose r = |x| and v 2 = v ■ v). 
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A. The multipole moments 



Our final result for the spin part of the mass-quadrupole moment at 2.5PN order (1PN 
order beyond the dominant SO term), for general orbits and in the CM frame, reads 
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(5.3) 



The scalar product of ordinary Euclidean vectors is indicated by parenthesis, e.g. (vS) = 
v • S, the cross product by the usual cross symbol, (x x X) 1 = e ijk x J Ti k , and the mixed 
product of three vectors by (S, x, v) — S • (x x v) = e l]k S l x 3 v k . We recall also that the STF 
projection is denoted using carets surrounding indices, i.e. (ij). Next, the spin part of the 
current quadrupole moment at 1.5PN order (also 1PN order beyond the leading term) is 
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(5.4) 



Notice that the 1.5PN current quadrupole moment sJij was also computed in Ref. [31], see 
Eq. (4.18) there. However our result (15. 4 ft differs from their result. There are two reasons for 
this discrepancy. The main reason is that Ref. 3J| completely neglected the non-compact 
support terms, which originate from the non-linearities of the Einstein field equations via 
the term in Eq. (12. ip . and physically represent the gravitational field acting as a source 
for the multipole moment. As we have seen in Sec. IIVI these terms are not negligible. Their 
contribution to the 1.5PN order current moment sJij has been computed in Eq. (|4.7bl) . The 
second reason for the difference between Eq. (15 ,4p and the result of [31] is a computational 
error in [31| when they apply the integration formula (I4.4p for computing the compact- 
support terms. We have already commented upon this error after Eq. (14.41) above. These 
two errors fully account for the discrepancy between our result and the one of Eq. (4.18) in 



Ref. 31 



Next, in order to derive the GW flux at 2.5PN order, we need the spin parts in the mass 
octupole and current octupole moments, but only at the lowest order in the spins. They can 
be obtained from our previous computation leading to Eqs. (13.61) . after the CM reduction. 
We find 6 



I ijk 
S 



J ijk 
S 



v |_9 Sm x (i xj , x s , k ) _ 3 (3 _ n v x<v(v x s)fe> 
{ 2 m 2 

+3 ^ zV(x x S f + 3 (1 - 3z/) x ( V(x x S) fe) | + O (±\ , 

-/2x<V5 fc >+2 — x^x^\+o(\) . 
c { m J V c / 



(5.5a) 
(5.5b) 



6 The result for s^ijk agrees with the one given by Eq. (4.17) in Ref. [31(- However, we notice a misprint in 



the second term of Eq. (4.17) in [31(, in which x^ k should read 
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B. The gravitational-wave energy flux 



With all these moments, Eqs. 05.3p -( ph4"]) and (I5.5p . and only with those, we can compute 
the 2.5PN spin part of the GW flux. Indeed, recall that the spins start at 1.5PN order 
in the mass moments and at 0.5PN order in the current ones, so one can easily see that 
in higher multipoles spins will enter the flux at higher PN order. On the other hand, one 
can check that it is not necessary to include the effects of tails of GWs, and more generally 
of any non-linear multipole interaction. Indeed, the tails give a correction to each of the 
source-type multipole moments II and Jl at the relative order 1.5PN ~ 1/c 3 (see e.g. jlsj]). 
For the mass quadrupole 1^ the spin itself is at order 1.5PN so the tail will arise only at 
order 3PN ~ 1/c 6 in the flux. For the current quadrupole Jy the spin is at 0.5PN but Jjj 
comes in the flux at 1PN order, so again we see that the corresponding tail will only be at 
3PN order in the flux. In conclusion, for this problem it is sufficient to express the flux solely 
in terms of the source multipole moments II and J^, all multipole interactions built in the 
radiative moments seen at infinity, namely Ul and Vl, are negligible. Furthermore, as we 
have seen only four multipolar contributions are important for this application. Therefore 
(cf. Eq. (4.28) in 0) 



T 



G 

75 



1 



I ij I ij + 



1 16 
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+ terms not contributing to the spins at 2.5PN order. 



(5.6) 



In order to compute time derivatives of the moments fl5.3l) - fl5.4p and ( 15. 5p . we must be 
careful at including the non-spin terms of and since these terms will generate by order 
reduction of the accelerations some new spin terms at 2.5PN order. In particular we need 
for this problem the non-spin part of the mass quadrupole 1^ at 1PN order, which is given 
by HEil, 
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The current quadrupole moment is also needed at 1PN order, 
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Note that for both 1^ and there are some contributions at 1.5PN order which depends on 
the spin variables and are generated from the Newtonian term evaluated in the CM; these 
contributions have been included in the results (I5.3l) - fl5.4p above. 

Finally, we obtain for the flux (in the general orbit case but in the CM frame) the structure 
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The non-spin pieces /n, /ipn, /i.spn and /2.5PN are already known and we shall need below 
the Newtonian and IPN terms which are given by [5(J 52| 
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Notice also that the non-spin terms /i.spn and /2.5PN include the contributions of GW tails. 
Here we do not deal with the spin-spin (SS) term at 2PN order which is given in Refs. 39l. 140 
We obtain the SO coupling part at 1.5PN order as 
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(5.11) 



and for this part we find perfect agreement with Kidder et al. [39|, |40j. Finally, for the 
next-order SO part our result is 
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VI. REDUCTION TO QUASI CIRCULAR ORBITS 

From now on we assume that when the binary enters the frequency bandwidth of the 
LIGO/Virgo/LISA detectors the orbit has been circularized by the gravitational radiation 
reaction effect. By circular orbit we mean an orbit which is circular when the gradual 
radiation reaction inspiral can be neglected, and when the effects of spins are averaged over 



time. With such proviso there is a well defined notion of a circular orbit (see Refs. |39l . 140 
and paper I). 

For circular orbits the orbital frequency u is linked to the distance r between particles 
in harmonic coordinates by a relativistic extension of Kepler's law, which has already been 
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given in paper I for what concerns the SO effects. Let us write it again here, but let us also 
add to it, for the benefit of potential users of these formulas, all the non-spin contributions 
up to the 2.5PN order, following known results from the literature ( e.g. J2l| and references 
therein). The 3PN and 3.5PN non-spin terms, computed in Refs. [20j, |23j, can be added 
straightforwardly if necessary. However, for convenience in this paper, we shall not display 
the non-linear spin-spin (SS) terms. Thus, all formulas of this Section will be complete up 
to 2.5PN order at linear order in the spins (i.e. but for the SS contributions). We have 
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(6.1) 



in which 7 = = O (c 2 ) denotes the harmonic-coordinate PN parameter. We recognize 
the lowest-order (1.5PN ~ 7 3 ^ 2 ) spin-orbit term and its 1PN correction at the 2.5PN ~ 7 5//2 
level. Here, as in paper I, we introduce an orthonormal triad {n, A, £} defined by n = x/r, 
£ = Ln/|Ln| where Ln = /xx x v denotes the Newtonian angular momentum, and A = 
£ x n. The quantities Sg and E^ in Eq. (16.11) are the components of the spin vectors (15. ip 
perpendicular to the orbital plane, namely Si = S • £ and E £ = S ■ £. The relation (16. ip 
can be inverted to give 7 in terms of an alternative PN parameter x, directly related to the 
orbital frequency through x = C^™^ ) 2 ^ 3 = O (c~ 2 ). As usual, it is better to express the PN 
formulas in terms of the frequency- dependent PN parameter x rather than 7 because they 
are invariant under a large class of gauge transformations. Hence, 
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The SO term at order 1.5PN ~ x ' is in agreement with Eq. (16) in [39 . 



The reduction of the GW flux J 7 , given by Eqs. (I5.11I) - (I5.12I) . to circular orbits is straight- 
forward, but care has to be taken from the fact that the non-spin parts of the flux at New- 
tonian and 1PN orders yield crucial contributions to the SO terms for circular orbits [beside 
the ones given by straightforward reduction of Eqs. (I5.11l) - (l5.12p ], Such contributions are 
generated by replacement of Eq. (16. ip into the 1PN flux given for general orbits by Eq. (I5.10p . 
Finally, we obtain 
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(6.3) 



or, equivalently, in terms of the PN parameter x, 
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All the non-spin terms are included up to 2.5PN order. Notice in particular the non-spin 
terms, proportional to 7T, which are at the same 1.5PN and 2.5PN orders as the SO effects; 

3/2 



these terms are due to GW tails |20l . |21J. For the leading SO term at order 1.5PN ~ x 
we find perfect agreement with Eq. (17b) in [391 ] . We also check that the result agrees in 
the test-mass limit zy — with the black-hole perturbation calculation of Tagoshi et al. [53| 
[see Eq. (G19) there]. 

The reduction of the center-of-mass energy E (computed in Sec. VII of paper I) to circular 
orbits is straightforward, and we simply report here the final result, completing it by the 
known non-spin terms 54|. We have 
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or, equivalently, 
E = 
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Alternatively, in terms of the single-spin variables the spin-dependent part of the above 
equation reads 
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where we denote by Ki = Si • £ for i = 1,2 the orientation of the spins with respect to the 
Newtonian angular momentum, and by Xi their magnitude defined in the standard way by 
Si = Gm 2 i Xi S { . 

Assuming non-precessing orbits, we list in Table [J the energy and the frequency at the 
so-called innermost circular orbit (ICO) (55j . The ICO is defined by the minimum of the 
center-of-mass energy for circular orbits expressed as a function of the orbital frequency u, 
and is computed from Eq. (16.61) . 
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TABLE I: Energy and angular frequency at the ICO for equal-mass [y = j) binary systems. The 
spins are maximal {xi = 1) an d have different orientations (k^ = 0, ±1). In three cases, indicated 
by a dash, there is no ICO, i.e. the energy function admits no real minimum. Spin-spin effects at 
2PN order are included. 



The orbital angular momentum (computed in paper I) in the case of circular orbits reads 
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For future use we give here the precessional equations evaluated in paper I, but reduced 
to circular orbits: 
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We recall our notation (15. ip for the spin variables. We denote by S n , S n and S\, T,\ the 
components of the spins along the vectors n and A respectively. As before we have neglected 
the SS terms. 

To compare easily with previous results in the literature, we have also computed the 
precessing equations for the spin variables Si and S 2 . They read 
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where S\\, Si n and S2A, ^Wi are the projections of the single-spin variables along A and n. 
We notice that with our choice of spin variables Si and S 2 the magnitude of the spin is 
not constant even when restricting Eqs. (16.111) to 1.5PN order. In Sec. IVHI we define some 
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alternative spin variables and S|, such that the magnitude of these spin vectors remains 



constant, i.e. • dS^/dt = with i = 1,2. The spin vectors agree with Kidder's (40 
spin variables at the IPN order, and generalize them to the next 2PN order. The main 
advantage of the definition is that the precession equations can be then written in the 
form dSi/dt — Oj x S^, where fij are the precession angular frequency vectors (given in 
Sec.lVill). 



VII. SPIN VARIABLES WITH CONSTANT MAGNITUDE 



In this paper and paper I we found convenient to use some specific spin variables Si 
and S2, defined in Sec. II of paper I. However, as discussed in paper I, other papers in 
the literature use a definition of the spin variables different from ours. For example, the 
spin-precession equations at IPN order in Ref. [40] read 
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where the superscript c stands for constant; in fact, the spin variables S^, S 2 are such that 
their norm or magnitude remains constant. Indeed one can readily check from Eqs. (17. ip 
that • dS1/dt = with % = 1,2. Our spin variables are related at the IPN order to the 
constant-spin ones (in the center-of-mass and for circular orbits) as 
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We can check that by taking the time derivative of the RHS of Eqs. f)7.2p . plugging in 
Eqs. (16. lip at IPN order, we recover Eqs. (17. ip . Note that the total angular momentum is 
invariant, since 

J = L + -Si + -S 2 = L c + -S; + -S c 2 . (7.3) 

c c c c 

Let us now define, at the 2PN order, in a general frame and for general orbits, some spin 
variables reducing to S\ and S 2 at the IPN order, and such that the magnitude of these 
spins remains constant. We shall still denote the latter 2PN spins as S^, S2; thus, we shall 
have, at the 2PN order, S£ ■ dS^/dt = with i — 1, 2. First of all we find that the new spin 
variables are related to the ones used in previous sections (and in the whole of paper I) by 
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together with the expression for S 2 obtained by exchanging all the particle's labels 1 -H- 2. 
In Eq. (17. 4p the notation is exactly the same as in paper I. The main advantage of such 
definition (17. 4p is that the precession equations can now be written into the form 



~dt 
dt 



Hi x , 
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(7.5a) 
(7.5b) 



showing that the spins precess around the directions of Oi and fl 2 , arid at the rates |f2i| 
and |f2 2 |. The precession angular frequency vectors fii, f2 2 can be computed up to the 2PN 
order by using the precession equations (6.1)-(6.3) of paper I. We find 
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together with 1 •f-)- 2. In the center-of-mass frame we obtain 
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For circular orbits, 
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Note that to obtain fi 2 we simply have to change 8m — > —5m. Recall that all these 
expressions are valid at linear order in the spins (excluding the SS terms); this means that 
with this approximation fix and fl 2 are independent of the spins. 

Finally let us express some of the main results of this paper in terms of the spin variables 
with constant magnitude. For the spin-dependent part of the circular-orbit energy [Eqs. (16. 6p 
or (16.71) ] we get 
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where S c e = S c • I and ^ = £ c • £, with S c = S c x + S| and £ c = - ^S c v Using 

the energy ( 17.9P as function of the constant spin variables, we have computed the ICO. 
For an equal-mass binary with spins anti-aligned with the orbital angular momentum, at 
2.5PN (3PN) order, we get E lco /m = -0.0122 {E lco /m = -0.0119) and mu lco = 0.064 
(muico — 0.061) to be compared with the numbers listed in Table HI The difference is not 
negligible. 

We also computed the spin-dependent part of the orbital angular momentum and the 
flux in terms of the spin variables with constant magnitude, 
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VIII. PHASE EVOLUTION AND ACCUMULATED NUMBER OF GW CYCLES 

In this Section we compute the time evolution of the binary's orbital frequency u, which 
results from the gravitational radiation reaction damping force. Instead of computing di- 
rectly the radiation reaction force, we use the standard energy balance argument 



T= rr 



dE 
dt 



U) 



where J 7 is the total emitted GW energy flux computed in Sec. |Vl and E denotes the binary's 
center-of-mass energy, namely the integral of the motion associated with the conservative 
part of the equations of motion, and which has been computed in paper I. Using F[x\ 
and E[x] expressed in terms of the spin variables with constant magnitude by Eqs. (I7.9P 
and (17. lip , we deduce Co from (18. ip . which is equivalent to 



to 

CO 
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2x 
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dx 



F[x] . 
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Notice that for this calculation it is crucial to use the variables associated with the constant 
magnitude spins, S c and £ c (rather than our original spin variables S and £ 8 ), since the 
spin variables S c and £ c are secularly constant, i.e. do not evolve by gravitational radiation 



We are grateful to A. Gopakumar and G. Schafer for pointing out this fact to us. 
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reaction. A proof that S c and S c are secularly constant (to the considered order) can be 
found in Ref. j56|. Hence the balance equation in the form of Eq. (18. 2p gives directly the 
secular evolution of the orbital frequency. 

During this computation the standard PN approximation is applied, i.e. we expand both 
the numerator and the denominator of Eq. (18.21) in the usual PN way, and finally express the 
result as a Taylor series in x. Other ways of addressing the computation, using particular 



PN resummation techniques, can be found in Refs. [14j , 1571 ] and references therein. We 
give the end result for the parameter £ = u/u 2 , which can be viewed as the dimensionless 
adiabatic parameter associated with the gradual inspiral, and which is dominantly of 2.5PN 
order (namely, the order of radiation reaction). In the final result, as everywhere else, the 
SO effects are at the 1.5PN and 2.5PN orders beyond the dominant approximation. We get 
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If necessary the non-spin contributions at orders 3PN and 3.5PN can be straightforwardly 
added 

Equations ([73]) . (EHD, (EHJ), f[7TTTj) and flOD with non-spin terms added through 3.5PN 
order and spin-spin terms included, together with the equation describing the rate of change 
of the orbital angular-momentum direction (deduced from L c = — -S5 — -^S^ at leading 
order) and the radiation field (see, e.g., Sec. C and Appendix B in Ref. |40j), can be solved 
semi-analytically, for special spin and mass configurations, or numerically. They provide 

19l and would need to be 



more accurate templates than currently used in the literature 
implemented for the search of GWs from spinning, precessing binaries. 

In the general case, taking into account the effect of precession of the orbital plane induced 
by spin modulations, the GW phase $gw is given by $ GW = G w + ^0gw, where 0gw is the 
"carrier phase" , defined by 0gw — 20 with <f> = f udt, and <50gw is a standard precessional 
correction, arising from the changing orientation of the orbital plane. The precessional 
correction <50gw can be computed by standard methods using numerical integration, see 
Ref. (9)]. Thus, the carrier phase G w constitutes the main theoretical output to be provided 
for the templates, and can directly be computed numerically from our main result, Eq. (18. 3p . 
In absence of orbital-plane's precession, e.g., for spins aligned or anti-aligned with the orbital 
angular momentum, the GW phase reduces to 0gw ; and the latter can be obtained by 
integrating analytically Eq. (I8.3p . We get 
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where 0o denotes some constant phase. In terms of the single-spin variables and S|, the 
spin-dependent part of the above equation reads 
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where xl an d ft? are defined by S£ = G mf xl and k\ = • The number of accumulated 
GW cycles between some minimal and maximal frequencies is 



GW 
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We list in Tables [TT] and II III the number of accumulated GW cycles (18. 6 p for typical binary 
masses in the most sensitive frequency band of ground-based and space-based detectors. For 
comparison we also show the contribution due to spin-spin terms at 2PN order evaluated in 
Ref. ' 
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20, 23 



as well as those due to the non-spin 3PN and 3.5PN orders computed in 
We denote £ c = S\ ■ S|. From the Tables HT1 and HTT1 we deduce two important results of this 
paper. First, we see that at 2.5PN order, if spins are maximal, i.e. \\ = 1, the number of 
GW cycles due to spin couplings is comparable to the number of GW cycles due to non-spin 
terms. Secondly, we find that for small mass-ratio binaries, the number of GW cycles due to 
linear spins at 2.5PN order can be much larger than the number of GW cycles due to spin- 
spin terms at 2PN order. These results thus show that the 2.5PN spin terms evaluated in 
the present paper have to be included in the GW templates if we want to extract accurately 
the binary parameters. 



TABLE II: Post-Newtonian contributions to the number of GW cycles (|8.6|) accumulated from 
<^min = vr x 10 Hz to w max = wisco = l/(6 3//2 m) for binaries detectable by LIGO and Virgo. For 
comparison, we add the contributions of spin-spin terms at 2PN order (we denote £ c = • S|) 
and non-spin terms at 3PN and 3.5PN orders. 





(10 + 1.4)M Q 


(10 + 1O)M 


(1.4 + 1.4)M 


Newtonian 


3577 


601 


16034 


1PN 


+213 


+59.3 


+441 


1.5PN 


-181 + 114^ + 11.8 4 xl 


-51.4 + 16.0 k\ x\ + 16.0 n\ x\ 


-211+65.7*5 Xl +65.7k§ xl 


2PN 


+9.8 -4.4 k\ n\ xl Xl + l-5 5 c x5 Xl 


+4.1 - 3.3 k\ K\ x\ X§ + 1.1 C c X\ X\ 


+9.9 - 8.0 k\ k% xl Xl + 2.8 ? xl Xl 


2.5PN 


-20 + 33.9 k\ x\ +2.9 k\ xl 


-7.1 + 5.7 kI xl +5.7k c 2 xl 


-11.7 + 9.3 nl xl +9.3 4 xl 


3PN 


+2.3 


+2.2 


+2.6 


3.5PN 


-1.8 


-0.8 


-0.9 



The number of accumulated GW cycles can be a useful diagnostic to understand the im- 
portance of spin effects, but taken alone it provides incomplete information. First, Agw is 
related only to the number of orbital cycles of the binary within the orbital plane, but it does 
not reflect the precession of the plane, which modulates the wave form in both amplitude 



and phase. These modulations are important effects. In fact, it has been shown [13|, U3, [16 
that neither the standard non-spinning-binary templates (which do not have built-in mod- 
ulations) nor the original Apostolatos templates [10| (which add only modulations to the 
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TABLE III: Post-Newtonian contributions to the number of GW cycles (|8.6p accumulated until 
w max = wisco = l/(6 3//2 m) over one year of integration, for binaries detectable by LISA. For 
comparison, we add the contribution of spin-spin terms at 2PN order (we denote £ c = S^ • S§) and 
non-spin terms at 3PN and 3.5PN orders. 





(10 6 + 10 6 )M© 


(10 6 + 10 5 )M© 


(10 5 + 1O 5 )M 


Newtonian 


2267 


4985 


9570 


IPN 


+134 


+281 


+323 


1.5PN 


-92.4 + 28.8 k\ x\ + 28.8 k\ x\ 


-243 + 161 k\ Xi + 11-5 nl xl 


-170 + 53 k\ x\ +53 4 X2 


2PN 


+6.0 - 4.8 k\ k\ xl xl + 1.7 ^ xl X% 


+12.5 - 4.4 k\ ^ Xi X\ + 1-5 ? c Xl X2 


+8.7 - 7.1 k\ k\ X \ X2 + 2-4 5 C Xi X2 


2.5PN 


-9.0 + 7.1 «5 Xi +7.1 K § x c 2 


-26.5 + 47.0k5 Xi +2.7k§ x\ 


-11.0 + 8.7«5xi+ 8.7 xl 


3PN 


+2.3 


+2.3 


+2.5 


3.5PN 


-0.9 


-2.3 


-0.9 



phase) can reproduce satisfactorily the detector response to the GWs emitted by precess- 
ing binaries. Modulations both in the phase and the amplitude of the wave form must be 
included 0, M, E3, M ■ Second, even if two signals have Aqw that differ by ~ 1, one can 



always shift their arrival times to obtain higher overlaps, but at the cost of introducing sys- 
tematic errors in the binary parameters. To quantify the impact of the 2.5PN spin terms in 
detecting GWs from spinning, precessing binaries, one should evaluate the maximized over- 
lap (fitting factor) between the 2.5PN template family and the 2PN template family used in 



Refs. [14|, |15|, |17|, uM- Those template families are defined by the GW signal computed along 
the binary evolution together with the spin and angular momentum precession equations. 
We expect that the maximized overlap between the 2.5PN and 2PN templates could be 
high, because the spins and the directional parameters entering the template families pro- 
vide much leeway to compensate for non-trivial variations of the phasing (see e.g., Table II 
in Ref. [l5| where maximized overlaps between several PN templates were computed). This 
study goes beyond the goal of this paper and will be tackled in future work. 



IX. CONCLUSIONS 

Within the multipole-moment formalism developed in Refs. (251-30]. we obtained the 
SO couplings, IPN order beyond the dominant effect, in the binary's mass and current 
quadrupole moments, as well as in the GW energy flux. The current-quadrupole moment 



with SO couplings at 1.5PN order was derived in Ref. [31], but our result differs from 
the expression computed there for two reasons, (i) The authors of Ref. [3l| neglected the 
non-compact support terms which originate from the non-linearities of the Einstein field 
equations and are not negligible at this order, (ii) Their result for the compact-support 
terms is affected by a computational error. The mass-quadrupole moment with SO couplings 
at 2.5PN order (including all compact and non-compact support terms) is computed here 
for the first time. 

The binary's energy and the spin precession equations including SO couplings through 
2.5PN order were computed in paper I. They were used to derive the secular evolution of 
the binary's orbital phase through 2.5PN order in the spins. We found that the 2.5PN terms 
give a relevant contribution to the number of accumulated GW cycles within the binary's 
orbital plane. In Tables [TT] and IHH we listed the number of GW cycles for typical binaries 
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detectable with ground-based and space-based detectors, such as LIGO/ Virgo and LISA. 
When spins are maximal, the SO contribution at 2.5PN order is comparable to that of 
the non-spin part at the same 2.5PN order. For some binary mass-configurations, the SO 
contribution at 2.5PN order can be larger than that of SS couplings at 2PN order. 

In order to extract accurately the parameters of maximally or mildly spinning binaries 
with ground-based detectors of first generation, having typical signal-to-noise ratio (SNR) 
of the order of 10, we expect the spin corrections through 3.5PN order to be sufficient. 
With space-based detectors having SNR of the order of 10 2 -10 3 , we would need a priori to 
compute non-spin and spin corrections at much higher PN order (for parameter estimation). 
For what concerns the impact of the SO couplings at 2.5PN order on the actual detection, 
it would be relevant to evaluate the maximized overlaps between templates that include SO 
effects through 2.5PN order against templates that include SO and SS effects through 2PN 
order. We anticipate that, at the cost of introducing systematic errors in the estimation 
of the binary parameters, the maximized overlaps could be high. In fact, the binary and 
directional parameters may compensate variations in the PN phasing. 

For future applications, we listed in Sec. IVIII the relevant equations defining the spinning 
dynamics and the GW phasing in terms of the constant spin variables. Such formulation 
is broadly used in the literature to define spinning, precessing templates for compact bina- 
ries @5a, H 



Finally, we computed the contributions of the spin terms to the location of the innermost 
circular orbit in the case of black-hole binaries. The results for equal-mass objects with 
maximal spins are summarized in Table HI Spin couplings at 1.5PN and 2.5PN orders can 
give significant contributions to the energy and frequency at the ICO (and nearby). 
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